Gravity-free hydraulic jumps and metal femtocups 
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Hydraulic jumps created by gravity are seen every day in the kitchen sink. We show that at small 
scales a circular hydraulic jump can be created in the absence of gravity, by surface tension. The 
theory is motivated by our experimental finding of a height discontinuity in spreading submicron 
molten metal droplets created by pulsed-laser ablation. By careful control of initial conditions, we 
show that this leads to solid femtolitre cups of gold, silver, copper, niobium and tin. 
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It has long been observed that water flowing hori- 
zontally can display a discontinuity in height [if 0, Q- 
This is the hydraulic jump, seen for example when wa- 
ter from a faucet impinges on the kitchen sink and 
spreads outwards. Gravity is a key ingredient in these 
well-understood large-scale hydraulic jumps, as discussed 
briefly below. In this Letter we show, remarkably, 
that the shallow water equations support solutions for 
a gravity-free hydraulic jump. The driver here is sur- 
face tension at the liquid-air (or liquid-vacuum) inter- 
face, and jumps may be expected to occur when relevant 
length scales are submicron. Our theoretical study was 
prompted by careful experiments showing that molten 
metal droplets impinging on a solid substrate display 
such a jump, solidifying into cup-shaped containers of 
femtolitre capacity. The droplets are created by laser- 
ablation of a solid metal target. Femtocups made of 
different metals on various substrates are formed under 
carefully maintained conditions of laser energy and sub- 
strate temperature. Outside this narrow range of pa- 
rameters we find what one would normally expect: the 
droplets solidify into lump-shaped structures on the sub- 
strate. The ability to make, and subsequently leach out, 
femtocups at will has potential applications ranging from 
nanoscale synthetic chemistry to single cell biology. 

Before describing the experiments and the femtocups, 
we discuss what causes a gravity- free hydraulic jump. 
Consider a steady axisymmetric jet of fluid of radius a 
impinging on a solid plate placed normal to the flow. 
The density of the surrounding medium is assumed to 
be negligible. The fluid then spreads radially outwards 
within a relatively thin film. The dynamics within the 
film is described by the axisymmetric shallow- water equa- 
tion 0,0, p 
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where r and z are the radial coordinate and the coor- 
dinate perpendicular to the solid wall respectively, with 
origin on the solid surface at the centre of the imping- 
ing jet. The respective velocity components are u and 
w. The total height h of the fluid above the surface is 



a function of r, and a prime thus denotes a derivative 
with respect to r. The parameters in the problem are 
the acceleration due to gravity, g, the surface tension co- 
efficient, (7 for the liquid-air or liquid vacuum interface, 
and the density p and the kinematic viscosity v of the 
impinging fluid. For incompressible axisymmetric flow 
the equation of continuity, in differential and in global 
form, reads 
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where Q = na 2 Uj is the steady inlet volumetric flow rate. 
A characteristic inlet jet velocity Uj is thus denned. It is 
reasonable to assume 3] a parabolic shape in z for the 
radial velocity 



u(r,z) = C(r)(z 2 - 2h(r)z) 
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satisfying the no-slip condition at the wall (z = 0) and 
the zero shear stress condition at the free surface (z = h). 
The analysis does not hinge on this assumption; any 
reasonable profile shape will give qualitatively the same 
results. Using Eq. (|2j) and the kinematic condition 
w =Dh/Dt = uh! at z = /i, the momentum equation (QJ 
integrated over z from to h reduces after some algebra 
to 
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where all lengths are scaled by a, and the 0(1) positive 
constant b = 2/5 for a parabolic profile. The left-hand 
side of (@J) represents inertia, and the three terms on the 
right hand side appear due to viscosity, gravity and sur- 
face tension respectively. The relative importance of the 
inertial term to each of these is quantified respectively 
by the Reynolds number R = uja/v R, the Froude num- 
ber F = u 2 /(ga), and the Weber number W = pu 2 a/a. 
In large-scale flows surface tension has been shown Q 
only to make a small correction to the location of the 
jump, so the last term is unimportant. This is to be ex- 
pected, since F in the kitchen sink is of order unity, while 
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W ~ 10 — 100. In contrast, consider Uj ~ 10 m/s and 
a ~ 10 _7 m, so F ~ 10 8 and W ~ 10 -2 . Here surface 
tension determines whether and where a jump will oc- 
cur, whereas it is the gravity term that may be dropped 
entirely from the equation. 

In general, a jump occurs if the pressure gradient be- 
comes increasingly adverse as the flow proceeds down- 
stream, and attains a magnitude large enough to counter 
the relevant inertial effects. The adverse pressure gra- 
dient may be created by gravity, or surface tension, or 
both. With gravity alone, Eq. Q reduces to 
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It is seen that if F is finite and h 3 r 2 is an increasing 
function of r, the denominator will go to zero at some r, 
i.e., a jump will occur in the framework of the shallow- 
water equations 5]. However, its precise location 
may not coincide with this estimate Q, and in a given 
experiment the radial extent available may be too small, 
or the inertia too low, for a jump to occur. Decreasing 
gravity has been shown to shift the jump location down- 
stream [l(J, consistent with Eq. (|SJ). Now considering 
surface tension alone, a crude prediction of the existence 
of a jump may be made by assuming the height upstream 
to be slowly varying in r, i.e., ti « 1, and thus setting 
h" = ti" = 0. We may then rewrite Eq. (gj as 
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A jump is now possible if h is an increasing function of r, 
which is a more stringent requirement than in the case of 
gravity. Note that the second term in the denominator 
appears due to radial spreading, i.e., surface tension alone 
cannot give rise to a one-dimensional jump like a tidal 
bore. 

We now solve Eq. as an initial value problem be- 
ginning at some location ri and marching downstream. A 
fourth-order Runge-Kutta algorithm is used. An initial 
radius somewhat larger than a is chosen, where it is as- 
sumed that a parabolic profile has been attained. The ini- 
tial conditions in h and its derivatives are not known ex- 
actly for this complicated problem, and numerical studies 
are being done to understand the flow in this vicinity. We 
have, however, repeated the computations with a variety 
of initial height profiles, and a range 1.2 < ri < 5 and 
0.1 < h < 1, and the results do not change qualitatively. 
Typical solutions are shown in figure |H At a particu- 



lar radial location r 



there is a singularity in the 



height of the fluid layer. Note that as we approach rj the 
shallow water equations are no longer valid, even approx- 
imately, so the present analysis cannot tell us anything 
about the actual shape close to or after the jump. The 
dependence of the jump location on the inlet jet veloc- 
ity uj is not monotonic, as seen from figure Et- Here 
the Reynolds number R is varied by changing Uj, with 
other quantities as in figure ^ so the Weber number in- 
creases as i? 2 , from 3 x 10 -9 to 90. For very low R or 
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FIG. 1: Typical solutions of Eq. 0), with F = oo, containing 
a singularity at r = rj. For demonstration, liquid properties 
are taken as those of molten silver (p = 5000 kg m -3 , v — 
10 -6 m/s 2 and a — 0.9Nm _1 ) and a — 5/xm. Solid line: Uj — 5 
cm/s; dashed line: Uj = 80 m/s. Values for molten tin show 
similar behavior. 
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FIG. 2: (a) The location rj of the singularity as a function of 
the inlet jet radius, expressed here in terms of the Reynolds 
number. For R = 0.01 - 90 (uj = 0.02 - 180m/s for the case 
considered) the jump radius is of the order of a few microns, 
as observed in the experiment described below, (b) Contour 
plot of jump location in the R — W plane. 



very high W , jumps are unlikely to form within the avail- 
able radius, i.e., inertia and surface tension must be in 
the right balance. The Reynolds and Weber numbers are 
now varied independent of each other (figure Eb). In the 
region shown in red rj > 60 so jumps are not predicted. 
(A higher cut-off does not change answers qualitatively.) 
Blue color indicates rj ~ ri, this region merits numeri- 
cal investigation. Gravity- free hydraulic jumps may be 
expected in the region shown by intermediate color, seen 
as a relatively narrow linear patch when R < 100. Here 
the jump location depends only on the ratio W/R. For 
a given W, jumps exist for over an order of magnitude 
variation in R. At R > 100, undular jumps are seen, 
which are being investigated further. 

We turn now to our experiments, which show a height 
discontinuity in spreading drops of molten metal. Since 
the experimental flow is transient in nature, a detailed 
comparison with the theoretical results is not possible, 
but the jump radius is in the right range. At larger 
scales hydraulic jumps are known to occur even when 
the incoming flow is in droplets rather than jets pj|. 
A Q-switched frequency tripled Nd:YAG laser (A = 355 
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FIG. 3: AFM image (right) and height profile (left) of a silver 
blob on silicon substrate kept at room temperature. 
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FIG. 4: Microscopy analysis of metal cups (a) SEM image of 
femtocups of silver on a silicon substrate obtained at E p = 
100 mJ/pulse and T s = 1173i^. A few blobs do exist, as 
indicated by arrows, (b) Typical height profile of a femtocup 
from AFM analysis, (c) Tilted field emission SEM image of 
a tin femtocup on silicon. 



nm, repetitive frequency, 10 Hz) is focused with pulse 
energy E p on a rotating metal disc in a vacuum cham- 
ber (10 -7 torr) and the resultant plume received at a 
distance of 4 cm on a clean vertical substrate held at 
a temperature T s , for a duration of 20 min 7]. The 
resulting metallic structures on the substrate are stud- 
ied by scanning electron microscopy (SEM), atomic force 
microscopy (AFM) and energy dispersive X-ray analysis 
(EDAX). Over most of the range of E p and T s , we ex- 
pect, and obtain, ill-shaped blobs of solidified metals, see 
figure El However, for a small range of these parameters, 
there is a strong preference to form cup-like structures 
of outer diameters ~ 300nm to 10/im, with side walls 
~ lOOnm high, and capacity ~ 1 fL (fig. QJ). The jump 
diameter is usually about half the total diameter. Height 
profiles associated with atomic force micrographs (figure 
as well as EDAX spectra (not shown) |jj confirm that 
the central region is raised from the substrate and con- 
tains metal. Interestingly, pulsed-laser ablation has been 
used extensively to produce a variety of structures 
but femtocups have not been reported before, although 
we notice stray instances of similar structures in other 
studies 0. 

This surprising femtocup structure is consistent with 
the proposed dynamics: of a droplet spreading out thinly 
initially and then undergoing a height discontinuity. We 
obtain femtocups of gold, silver, copper, tin and nio- 
bium of repeatable statistics on glass, silicon and graphite 



FIG. 5: SEM images with different metals, laser energy and 
substrate temperature, (a) Cu femtocups on Si. E p =100 
mJ/pulse, T s = 300 K; (b) Nb on Si, £ p =100 mJ/pulse, 
T s = 1273 K; (c) Ag on Si, E p = 60 mJ/pulse. The number 
density of well-formed cups is lower, (d) Ag on Si, E p = 100 
mJ/pulse, T s = 773 K much less than T m = 1234 K. The cups 
are not well- formed. Inset, T s = 1273 K > T m . Only patches 
are observed. 



(HOPG), see examples in figure 0i and b, and later in 
figure [7| The solid surface being vertical and the length 
scales small mean that the effect of gravity is negligible. 
Inertia on the other hand is considerable, since velocities 
are high. We do not have a direct estimate of Uj, but we 
may estimate it from earlier measurements in many sim- 
ilar experiments to range between lm/stolOOm/s. 
Also, the range of R and W over which a surface-tension 
driven hydraulic jump occurs translates to a particular 
range of E pi since laser fluence determines scales and 
speeds in the incoming jet, compare figures 0Ji and 03. 
While the substrate is hotter than the metal's melting 
point T m , the cup may form initially but cannot solidify 
and liquid flows back into the cup, so the final object 
is as seen in the inset of figure 01. With T s far below 
T m the tendency to form cups is much reduced (figure 
01, probably because solidification is too rapid for flow 
to be completed. Optimal conditions are thus E p (~ 100 
mJ/pulse for silver) and T s close to but below T m . Out- 
side the correct range, blobs form rather than cups. 

That the jump is directly related to droplet dynamics 
is confirmed by varying the substrate orientation with 
respect to the incoming jet (see schematic in Fig. |BJ. As 

the inclination of the substrate away from the normal, 
is increased, the structures become increasingly elliptical, 
especially beyond 40°, in accordance with the azimuthal 
variation of R and W. 

Since the experiment includes additional complexity in 
the form of solidification, we estimate relative time-scales 
of jump formation tj and solidification of a droplet t c . For 
the experimental values of substrate thickness, t c ranges 
from ~3x 10 _4 s on silicon to ~ 10 _2 s on glass (taking 
into consideration conduction, radiation and latent heat), 
while tj ~ rj/uj ~ 10 _6 s or less. Contact-line freez- 
ing can give rise to an increase in height in the vicinity, 
typically amounting to a small percentage of the height 
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FIG. 6: Elliptical femtocups with inclined jets. Inset: Sample 
SEM image of tin obtained at = 40 ± 1.5°, denning a and b. 
The ratio b/a increases with 0, each data point is an average 
on several cups. The major axis lies along the direction of 
maximum flow (dashed arrows). 
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FIG. 7: Histograms and SEM images of tin femtocups on (a) 
glass and (b) silicon, deposited simultaneously. 



at the contact line [ljj. We do notice a dependence on 
the substrate of the size distribution of femtocups (fig- 
ure 0) and also some visual differences in the shape of 
the femtocup. Our ongoing numerical study, including a 
non- uniform temperature profile and its effects, is there- 
fore aimed at a better representation of the experiments. 
Also being addressed are the experimental finding of ra- 
dial striations in the femtocups under certain conditions, 
and the theoretical finding of undular hydraulic jumps 
(similar to in other conditions) when R is greater 
than about 25/ W (figure®. 

In summary, for the first time, a hydraulic jump solely 
driven by surface tension is shown to occur. Experimen- 
tally we show evidence for such jumps in submicron high 
inertia droplets of molten metals spreading radially out- 
wards on a substrate. The detailed shape in the vicinity 
of the jump and the transient problem including the so- 
lidification process is being studied numerically. 

We are grateful to Prof. G Homsy, Prof. CNR Rao 
and Prof. R Narasimha for useful discussions. RG and 




FIG. 8: (a) Striations in the structure, (b) An undular hy- 
draulic jump at R = 5000, W = 0.005. 



in the central region. Contrast this to our jump where 
the height at the rim is several-fold larger than that in 
the central region. In spite of this, and the disparity in 
time scales, we cannot rule out a role for local freezing 
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